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Abstract. A variety of descent and major-index statistics have been defined for symmetric groups, 
hyperoctahedral groups, and tfieir generalizations. Typically associated to pairs of such statistics is 
an Euler-Mahonian distribution, a bivariate generating function identity encoding these statistics. 
We use techniques from polyhedral geometry to establish new multivariate generalizations for many 
of the known Euler-Mahonian distributions. The original bivariate distributions are then straight- 
forward specializations of these multivariate identities. A consequence of these new techniques are 
bijective proofs of the equivalence of the bivariate distributions for various pairs of statistics. 
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1. Introduction 

The symmetric group Sn is the group of permutations of [ra] := {1, 2, . . . , n}, also realized as the 
Coxeter group ^n-i yielding symmetries of a simplex. For a permutation tt & Sn, the descent set 
is a classical object of study in combinatorics. 

Definition 1.1. Let tt G Sn- The descent set of tt is 

Des(7r) := {j G [n - 1] : 7r(i) > 7r(j + 1)} . 

The descent statistic is des(7r) := #Dcs(7r). 

The descent statistic is encoded in the Eulerian polynomial X^TreSn t'^'^^^'^^ and the most basic 
identity for Eulerian polynomials is 

^des(7r) 

A;>0 '''' 

Euler used this identity to define Eulerian polynomials in [13] which he needed in his study of what 
is now called the Riemann (^-function; one can safely assume that he was aware of the connection 
of his polynomials to descent statistics (see [15] which also gives connections between Eulerian 
polynomials and algebraic geometry). 

Equation (1) has inspired a host of generalizations and extensions. The first such extension is the 
following ^'-analogue of (1), often attributed to Carlitz [10] but actually going back to MacMahon 
[19, Volume 2, Chapter IV, §462]. This extension involves a joint distribution of the descent statistic 
and the major index, defined as follows. 

Definition 1.2. For tt G Sn, the major index of tt is 

maj(7r) := ^ j. 

jGDes(7r) 

MacMahon's extension of (1) is the following, where [m]q := 1 + q + + ■ ■ ■ + q"^~^. 
Theorem 1.3 (MacMahon). 



r, ides(7r)„maj(7r) 



k>0 

Note that (1) follows from Theorem 1.3 by setting q = 1. The search for further generalizations of 
this identity, called the Euler-Mahonian distribution or the Carlitz identity, has focused on finding 
new identities of the form 

^stati(g) stat2(g) 

(2) E[fe + 1]^'= hitq) 

k>o LLj=o'T-j[hq) 

for various families of groups G„ and statistics stati and stat2 defined on elements of Gn, together 
with naturally occuring families of functions hj{t, q). This search has been successful, also produc- 
ing analogous generalizations of the identities (6) and (7) discussed in the next section. To our 
knowledge, there are three general approaches to proving such identities: 

• via combinatorial/bijectivc proofs in the theory of partitions and their extensions; 

• via connections between permutation statistics and the theory of Coxeter groups, including 
connections to invariant theory and the coinvariant algebra of a Weyl group; and 

• via the theory of symmetric/quasisymmetric functions. 
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For more information regarding the first two approaches, see the citations hstcd throughout this 
paper. For examples of the symmetric/quasisymmetric function approach, sec [17, 20, 23]. 

Our goal is to provide new multivariate generalizations of these identities using polyhedral geome- 
try and lattice-point enumeration; as a consequence, we obtain new proofs of two- variable identities 
in the form of (2). One of the benefits of the geometric approach is that it is relatively simple, the 
key ingredients being the triangulation of the unit cube by the braid arrangement together with 
careful choices of ray generators for unimodular cones. Another benefit is that bijcctivc proofs of 
the equidistribution of various pairs of statistics are obtained as immediate corollaries. 

As we discuss in Remark 4.3, our multivariate identities can be viewed as Hilbert series identities 
for various finely-graded algebras, i.e., algebras equipped with an N"-grading. A Hilbert series ap- 
proach to multi-variate extensions of these identities has previously been used in [2] and subsequent 
papers, emphasizing the use of descent bases for coinvariant algebras. Our algebras and special- 
izations are in some sense more straightforward than the previously considered ones, because the 
bivariate identites arise as simple specializations of our multivariate identites, requiring minimal or 
no additional substitutions and algebraic manipulations. The geometric perspective also allows us 
to avoid the use of straightening laws and other algebraic techniques regarding coinvariant algebras. 

Our paper is structured as follows. In Section 2, we discuss analogues of (1) for generalizations of 
permutation groups. In Section 3, we discuss the results we will need from integer-point enumeration 
and polyhedral geometry. In Section 4, we use polyhedral geometry to prove Theorem 4.1; this 
proof serves as a model for all the proofs in the paper. We also briefly discuss connections between 
our approach, the theory of P-partitions, and the theory of affine semigroup algebras. 

Sections 5, 6, and 7 contain our new results. In Section 5, we prove Theorems 5.7, 5.9, 5.14, 
and 5.15. These four theorems give new multivariate identities for negative and flag statistics in 
type B. In Section 6, we prove Theorems 6.3 and 6.4, new type-D generating function identities. 
Finally, Section 7 contains theorems extending our type-B results to wreath products of cyclic 
groups with symmetric groups. 

2. Generalized permutation groups and descents 

We discuss in this section analogues of (1) for hyperoctahedral groups, type-D permutation 
groups, and wreath products of cyclic groups with symmetric groups. 

2.1. Hyperoctahedral groups. The hyperoctahedral group Bn, the Coxeter group yielding sym- 
metries of a cube, consists of signed permutations [22], i.e., pairs (tt, e) where tt E Sn and e G {±1}". 
For elements of there are several deflnitions of descents in the literature; we provide three of 
them here. 

Definition 2.1. For an element (tt, e) G -B„, the naturally ordered descent set is 

(3) NatDes(7r, e) := (j G {0, 1, . . . , n - 1} : ej7T{j) > ej+nr{j + 1)} , 

with the convention eo7r(0) = 0. The natural descent statistic for Bn is natdes(7r, e) := #NatDes(7r, e). 

The reason for calling this the naturally ordered descent set is that it uses the natural order 
— n < —n -|-l<---<— l<l<2<---<non the integers. Thus, the permutation (written as a 
word) 

-3 2-1-4 

in B4 has descents in the second and third positions. The next two deflnitions of descent set for i?„ 
are special cases of descents in the context of colored, or indexed, permutations arising from wreath 
products. 
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2.2. Wreath products. The wreath product I Sn of a. cychc group of order r with Sn consists 
of pairs (tt, e) where tt G S'n and e G {u}^,u}^, . . . ,uj^~^}^ for lo := e^'^*/''. Thus, e is a sequence of 
powers of an r^^ root of unity; the group Bn arises as the wreath product Z2 I Sn- By convention, 
for elements of Z,. I Sn we define additional values of tt and e as follows: iTn+i := n + 1, e„+i := 1, 
TTo := 0, and eo := 1. 

We will find it convenient to use window notation for elements of wreath products: If €j = lo'^^ , 
then we will denote (vr, e) as the window [7r(l)'^^ 7r(2)'^2 ... 7r(n)'^"]. Wc use the notation j^J and 
(w^J , j) interchangeably for elements of {uj^,uj^, . . . ,u^~^} x [n]. It is sometimes convenient to refer 
to -7T{jY^ as 7r(j) with color cj. In [27], Steingn'msson defined the following descent set for elements 
of 7dr I Sn- 

Definition 2.2. Totally order the elements of {oj^ , lo^ , . . . , lu'^~^} x [n] by j"^^ < k'^'' if Cj < Ck 
or if both cj = and j < k hold. For an element (tt, e) = [7r(l)"^i 7r(2)"^2 . . . 7]-(7^)c„j ^ 
Steingrimsson's descent set is 

(4) StDes(7r,e) := {j € {1, . . . ,n} : TT{jf^ > 7r(j + lY^+^]- 
Steingrimsson's descent statistic is stdes(7r, e) := #StDes(7r, e). 

As an example, observe that with Steingrimsson's ordering we have {uP ,uj^ jUp'} x [3] ordered as 

1° < 2° < 3° < 1^ < 2^ < 3^ < 1^ < 2^ < 3^ 

and the permutation 

[223^1!] 

has descents in positions 2 and 3. 

Finally, we define the following closely-related descent set. This definition differs from Ste- 
ingrimsson's both in the role played by the order of the roots of unity and in the indices where 
descents may occur. 

Definition 2.3. Totally order the elements of {uj'^~^ ^uf'"^ , . . . ,0;*^} x [n] by < k'^'^ if Cj > Ck or 
if both Cj = Ck and j < k hold. For an element (tt, e) = [7r(l)'^i Tr{2y^ - - - 7r(n)"^"] in I Sn, the 
descent set is 

(5) Des(^, e) - {j e {0, . . . , n - 1} : tt^T^ > 7r(j + 1)^^+^ }. 
The descent statistic is des(7r, e) := #Des(7r, e). 

As an example, observe that with this order we have {a;°,a;^,a;^} x [3] ordered as 

1^ < 2^ < 3^ < 1^ < 2^ < 3^ < 1° < 2° < 3° , 

and the permutation 

[32 2° li] 

has descents in positions and 2. 

The Eulerian polynomials for wreath products are given by e)eiriS„ t'^^^^'^''^^ where one may 
use either of the two wreath product descent definitions or, in the case r = 2, the natural descent 
statistic. The resulting analogue of (1) is 



(6) ^(rA; + lfi 



Ej.des(7r,e) 
(7r,e)6ZwS„ ^ 



(1 - 

fe>0 

This identity appears to have been found by various authors for different descent statistics; for 
more details, see [9, 27]. 
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2.3. Type D. One may also consider type--D Eulerian polynomials stemming from the signed 
permutations in with an even number of — I's. Let 

Dn ■■= {(tt, e) € Bn'. ei • • • e„ = 1} . 

The definition of DNatDes(7r, e) and dnatdes(7r, e) in type D is analogous to (3), except that we 
now use the convention eo7r(0) := — e27r(2). Brenti [9, Theorem 4.10] proved that 

^dnatdes(7r,e) 

(7) ^ {{2k + ir - 2"-^ {Bn{k + 1) - Bn{0))) = ^"'^|f:,.„+i , 

fc>0 ^ ^> 

where Bn(x) is the n'th Bernoulli polynomial. 

3. A GEOMETRIC PERSPECTIVE 

The forms of equations (1), (6), and (7) suggest that one should look at them geometrically as 
stemming from lattice-point enumeration of the cube [0, r]" as it is partitioned in various ways; 
for example, (1) suggests we consider [0, 1]" partitioned by the braid arrangement consisting of the 
hyperplanes Xj = Xk ior 1 < j < k < n. We will encounter certain simplices below, which generally 
look like (after a change of variables) 

Ai:=LeR^: ^ < < x^ ,<■■■< x, < 1, 

[ xj+i <xj it J e I 

where / C [n] is some index set, and we use the convention > if n G /. The vertices of A/ are 

0, ei H h Sn, ei H + e„_i, . . . , ei + 62, ei, where ej is the j'th unit vector in R". Note that 

A/ is unimodular, i.e., the n edge directions at any vertex of A/ generate Z". 

The cone over A/ is the nonnegative span of {(l,a;) G M""*"^ : x G A/|, where we encode the 
"new" dimension by the variable xq, i.e., 

cone (A/) := R>o cq + ^ M>o (eo + ei + 62 H h ej) + ^ M>o (cq + ei + 62 -I h ej) , 

jei j^i 

where the complement of I is taken in [n] . Let 

ac{zo,zi,...,Zn) := ^ z"^ 
meCnZ"+i 

be the multivariate ("full") generating function encoding the integer lattice points in a cone C, 
where we have used the shorthand z"^ := zj^^z^^ ■ ■ ■ z^". A standard geometric-series argument 
(see, e.g., [5, Chapter 3]), together with the unimodularity of cone (A/), gives the following. 

Lemma 3.1. Let A/ be as above. Then 



^cone(/\T)\ZQi Zl, . . . , Znj — \ 

\\j=0V-- Z0ZlZ2---Zj) 



li=0 ' 

Wc will not always use the above natural way to write the generating function of a unimodular 
cone, in which case we will apply the following more general lemma (for details, again see, e.g., [5, 
Chapter 3] or [26, Section 4.6]). 

Lemma 3.2. Let C = X]j=o^>o'*^j + Z^j=fc+i ^>oi'j a half-open simplicial cone in M."'^^ with 
linearly independent generators Vq, Vi, . . . ,Vn G Z""^-*^. Then 

an{zo,zi,...,Zn) 



ac{zo,zi,...,Zn) 



n -=0(1-^^0 
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where li := Ei=o[0' + E}=fe+i(0, l]vj. Furtherm, ore, the number of integer points inU ( and 
thus the number of monomials in auizo, zi, . . . , Zn) ) is given by the determinant of the matrix with 
column vectors Vq, Vi, . . . ,Vn- 

We refer to 11 as the fundamental parallelepiped of C; note that it depends on the choice of 
generators of C. 

4. Type A 

We begin with a multivariate identity that specializes to Theorem 1.3. The proof of this identity, 
though simple, demonstrates the approach used in this paper. 

Theorem 4.1. 

TTfyt + 1]^ Z^= y nj£DesW^0^^(l)M2)---^^(j) 

^0 ° n;=o (1-^0 ^.(1) ^.(2) • • • z^j)) ■ 

Proof. Triangulate the n-cube [0, 1]" into the disjoint union of simplices 

A. := L G M" : ° " ^^^^^ " ^-(""'^ "t^ ' ' A"""^'^ ' ^' 
i ^7r(i+i) < if J G Des(7r) 

(one for each vr G S'n)- The strict inequalities determined by the descent set of tt ensures that this 
triangulation is disjoint. By Lemma 3.1, 

. ^ _ njsDcs(7r) ^0 ^7r(l)^7r(2) ' ' ' ^7r(i) 

<''cone(A^)l^05 2:1, ... , Zn) — / x • 

llj=0 (1 - ^0 ^;7r(l)^;7r(2) • • • ^7r(i)j 



On the other hand, 



O-cone([0,l]") (^^0, Zl, . . . , Zn) = (l + Zj + z] -\ h Zj^ Zq , 

fc>0 j=l 

and the disjoint triangulation gives 

'''cone ([0,1]") {zo,Zi,...,Zn) = ^ 

'-'^cone(A,r 

){zo, Zi, . . . , Zn) . □ 

neSn 

Proof of Theorem 1.3. Setting t := zq and q := zi = Z2 = ■ ■ ■ = Zn va. Theorem 4.1 gives 

Remark 4.2. Our approach is related to the theory of P-partitions [24, 26], which can be interpreted 
geometrically through triangulations of the unit cube. Given a finite poset P, the set of P-partitions 
can be interepreted as the integer points inside a cone admitting a nice unimodular triangulation, 
where the unimodular subcones are indexed by linear extensions of P. Our approach is based 
almost entirely on the triangulation of [0, 1]" induced by the braid arrangement; the relationship 
with P-partitions is that [0, 1]" is a truncation of the P-partition cone in the case where P is an 
antichain of size n. Thus, [0, 1]" is the order polytope of the antichain of size n. As a result, our 
Theorem 4.1 is similar to [24, Theorem 7.1]. However, in the rest of the paper, we repeatedly choose 
non-unimodular generators for unimodular cones to obtain our generating function identities. To 
our knowledge, this approach has not been used to study P-partitions. 
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Remark 4.3. The generating function in Theorem 4.1 is the finely-graded Hilbert series for the affine 
semigroup algebra with semigroup given by the set of all integer points in the cone over [0, 1]", as 
discussed in [14, 21, 25]. The study of semigroup algebras arising from polyhedral cones has been 
an area of intense study for combinatorial commutative algebraists over the past several decades. 
The most important general result regarding Hilbert series for such cones is Hochster's theorem, 
which states that normal affine semigroup algebras are Cohen-Macaulay [16] . The Cohen-Macaulay 
property forces serious constraints on single-variable specializations of the associated finely-graded 
Hilbert series for the algebra. 

5. Type B 

5.1. g-analogues involving (k + 1)"^. In the case Gn = -B„, there are several pairs of statistics 
resulting in (j'-ifications of (1). In [1], Adin, Brcnti. and Roichman introduced a variety of relevant 
statistics. We first recall a descent statistic for elements of that reads a word as if it were a 
permutation in the symmetric group. 

Definition 5.1. Use the order — 1 < • • • < —n < 1 < • • • < n on [— n,n] \ {0}. For an element 
(tt, e) G Bn we define the type-A descent set as 

DesA(7r,e) := {i e [n - 1] : eiiTi > ej+iTTj+i} 

and the type-A descent statistic as 

desA(7r, e) := #DesA(7r, e) . 

The type-A major index is 

major^(7r,e) := ^ i. 

jgDesx(7r,e) 

The second set of statistics we need are referred to as the negative statistics for type B. 

Definition 5.2. For an element (vr,e) S Bn, we define the negative set of (tt, e) to be 

Neg(7r,e) := {i G [n] : = -1} 

and we define neg(7r, e) := #Neg(7r, e). For a sign vector e, we denote by neg(e) the number of — I's 
in e. We define the negative descent multiset as 

NDes(7r, e) := DesA(7r, e) \J{nii) : = -1} = DesA(7r, e) [j Neg((7r, e)-^) . 

The negative descent statistic is 

ndes(7r, e) := T^NDes(7r, e) 

and the negative major index is 

nmaj(7r,e):= ^ i. 

«eNDes(7r,e) 

Note that when summing to obtain the negative major index, one sums over a multiset, so there 

is possible repetition of the summands. Originally, the negative statistic was defined using the 
natural order, but the following theorem holds for the negative statistics defined with either order. 

Theorem 5.3 (Adin-Brenti-Roichman). 

Ej.ndcs(7r,e) „nmaj(7r,e) 

h (i-on-=i(i-iVo ■ 

The final set of statistics obtained by Adin, Brenti, and Roichman are called the flag statistics. 
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Definition 5.4. Use the order —1 < • • • < — n < 1 < • • • < n on [— n,ra] \ {0}. For an element 
(tt, e) G Bn, define the flag descent statistic to be 

fdes(7r, e) := 2 • desyi(7r, e) + ci , 

where ci G {0, 1} satisfies ei = (—1)^^. The flag major index is 

fmajor(7r, e) := 2 • major^(7r, e) + ncg(7r, e) . 

The flag statistics were originally presented using the natural order, in which context the flag 
major index is denoted by fmaj rather than fmajor. As discussed in the remark in [1, p. 218], the 
following theorem holds for the flag statistics defined with either order. 

Theorem 5.5 (Adin-Brenti-Roichman). 

EE j.fdcs(7r,e) fmajor(7r,e) 
It. _L 1 1'-f"- — (^.e)6-B„ y 
,J " (i-on^=i(i-*v^) ■ 

5.2. Multivariate identities. We state and prove in this subsection multivariate generalizations 
of Theorems 5.3 and 5.5. Our multivariate generalizations are different than those obtained in [2] 
via coinvariant algebra theory. 

Definition 5.6. Define the subset of increasing elements of Bn, denoted l2,n, to be those elements 
satisfying desA{p, e) = 0, i.e., 72,n contains all permutations (p, e) satisfying the following conditions: 

• Cj = — 1 for J < neg(e) and ej = 1 for j > neg(e); 

• < P(2) < ••• < p(neg(e)); 

• p(neg(e) + !)<•••< p(n). 

It is straightforward that every element of Bn can be represented uniquely as (p, e)7r — (p^r, e7r) 
for some ir & Sn and (p, e) G l2,n, where pn is the standard composition of permutations and 

Ctt = (e7r(l)) e7r(2), ■ ■ ■ ^77{n))- ThuS, 

Bn= [J h,nT^ ■ 

Finally, observe that 

(8) Neg([(p,6)7r]-i) = Neg((p,e)-i) = {p(i) : ej = -1}. 

Theorem 5.7. 



k>Qj=l 



y-v j€Des(7r) j6Neg((p,e)-i) 



7re5„ (p,6)e72,„ (1 _ zq) n - z^zz,,-, ■ ■ ■ z. 



2 2 2 
■0^7r(l) ■■■^7r(i) 



Proof. This proof proceeds in two stages. We first triangulate the cube [0, 1]" into a disjoint union of 
simplices, then set up an indexing system for the integer points in the fundamental parallelepipeds 
for the cones over these simplices. Second, we bijectively associate the elements of Bn with these 
integer points in a way that allows us to recover, in our subsequent proofs, the descent and flag 
major statistics from these integer points. 
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We begin with the triangulation of cone([0, 1]") into the set of cones {cone(A7r) : tt G Sn} found 
in the proof of Theorem 4.1. While cone(A7r) is unimodular for each tt, for this proof we use the 
non-unimodular ray generators 

eo, 2(eo + e^(i)), 2(eo + e^(i) + e^(2)), . . . , 2(eo + e^(i) H h e^(„)) . 

There are 2" integer points in the fundamental parallelepiped for cone(A7r) using these ray gen- 
erators. Each such point can be expressed as a linear combination of the last n generators with 
coefficients aj € {0, plus a sum of shifting vectors for those integer points that need to be shifted 
away from the boundary of the cone. Typically one would use ray generators as shifting vectors 
for only those vectors that need to be shifted off the boundary. Due to the fact that cone(A7r) is 
unimodular, in this case we may instead shift all the points in the fundamental parallelepiped by 
the minimal ray generators bq + 6,^(1) + • • • + when j S Des(7r). This can also be justified by 
noting that the fundamental parallelepiped for cone(A^) with our non-unimodular ray generators 
is lattice-isomorphic to a half-open cube of the form [0, 1) x [0, 2)"'. 

Thus, every point p in the (shifted) fundamental parallelepiped for cone(A7r) can be uniquely 
expressed as 

n 

j6Des(7r) j=l 

with aj € {0, ^}. Associate to the point p the element (p, ejn G Bn, where aj = ^ if and only if 
j G Neg((/3, e)~"^) = Ncg([(p, e)7r]^^). This correspondence creates a bijection between the elements 
of Bn and the (appropriately shifted) integer points in the fundamental parallelepipeds for the 
cones over the A^. Note that this bijection encodes /2,n as the integer points in the fundamental 
parallelepiped for cone(Ai(j), where Id denotes the identity element in Sn- 
Thus, we have that 

JJ ^0Zn{l)Zn{2) ■ • ■ ^ttQ-) JJ ^0^7r(l) ^7r(2) ' ' ' ^7r(i) 
j6Des(7r) jeNeg((p,e)-i) 
0"cone(A^)l^O, ■ ■ ■ , Zn) - 2^ H • 

M^^^^r. (i-^o)n(i--o4(i)---4o)) 

This completes our proof, since from our triangulation it follows that 

O'cone([0,l]")(^0, ■ ■ ■ , Zn) = ^ ){zo,... 

1 Zn) ■ I— ' 

7re5n 

Proof of Theorem 5.3. Setting t = zq and q = zi = ■ ■ ■ = Zn m Theorem 5.7 yields our desired 
form on the left-hand side of our identity, while the denominator of the right-hand side uniformly 
becomes 



(i-t)n(i-iv^)- 



Each element (p, e)7r G ^2,n7r contributes to the numerator on the right-hand side of our identity 
a summand of 

n n 

jeDes(7r) jeNeg([(p,e)7r]-i) 
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Because of (8) and the fact that Des(7r) = DeSyi((p, e)7r), it follows that 

n ^^'^ n tgJ' = t'i<ies(p(7r,€))^nmaj(p(7r,€))^ 

jeDes(7r) jeNeg([(p,€)7r]-l) 

hence our proof is complete. □ 
Definition 5.8. For an element e G {±1}'*, set e„+i := 1 and for j G [n] set 

^ ■ \ else. 

Define ch(s) := G [n] : aj = 1}, the number of sign changes in e. 
Theorem 5.9. 

n 

n ^0^1(1)^1(2) ■ ■ ■ ^l(j) H Z0ZTr(l)Zn(2) ' ' ' ZTr(j) 

Proof. As in our previous proof, this proof proceeds in two stages. The key idea in this proof, as 
well as in the proof of Theorem 5.5 that follows, is that the encoding of an element (tt, e) G as 
an integer point is done using not the signs of e, but the sign in e given by the a^'s. 

We begin with the triangulation of cone([0, l]") into the set of cones {cone(A7r) : vr G Sn) found in 
the proof of Theorem 4.1. As in the proof of Theorem 5.7, for cone(A7r) we use the non-unimodular 
ray generators 

eo, 2(eo + e^(i)), 2(eo + e^(i) + e^(2)), • • • , 2(eo + e^(i) H h e^(„)) . 

In this case, we do not want to shift every integer point in the fundamental parallelepiped for 
cone(A7r) by a uniform shifting vector; we instead only want to shift the points on the boundary of 
the cone, and those by the appropriate non-unimodular ray generators. Hence we represent every 
such point p uniquely with a vector a G {0, of coefficients for 

n 

p = E 2aj(eo + e^(i) + • • • + e^(j)) + ^ 2(eo + e^(i) + • • • + e^(j)) . 

j=l jeDes(7r) 

a,=0 

We may then associate to the point p the element (vr, e) G B^ where vr is the same as the index on 
cone(A7r) and e is defined by = 2a j. This bijectively relates Bn to the (possibly shifted) integer 
points in the fundamental parallelepipeds of the cones over the A^^'s. Thus, we have that 



e(A^){zQ,...,Zn)= E 



jeDes(7r) 3-0'%=^ 



c^'^)^^" (i--o)n(i--o4(i)---4o)) 
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(Note that in the sum on the right-hand side, tt is fixed while e varies.) This completes our proof, 
since from our triangulation it follows that 

'7cone([0,l]")(^0, ■ ■ ■ , Zn) = ^ t^cone(A^) (^0, ■ ■ ■ , Zn) ■ □ 

7re5„ 

Proof of Theorem 5.5. Setting t = zo and q = zi = ■ ■ ■ = Zn in Theorem 5.7 yields our desired 
form on the left-hand side of our identity, while the denominator of the right-hand side uniformly 

becomes 

n 

(1 - 1) n (1 - tv^) ■ 

Each element (tt, e) G S„ contributes to the numerator on the right-hand side of our identity a 
summand of 

n i'l'' n 

jeDes{n) j;a^=l 

Therefore, our proof will be complete once we prove that 

(9) fdes(7r,e)= ^ 2+ E ^ 

jeDes(7r) j-a)=l 
a5=0 

and 

(10) fmajor(7r, e) = ^ 2j + ^ • 

jeDes(7r) j:a^=l 

A descent in position j of (tt, e) can arise for one of three reasons: either ej = 1 while e^+i = — 1, 
or ej = ej+i and j G Des(7r), or j = 1 and ei = —1. When viewing (tt, e) in window notation and 
reading right to left, every descent of the first kind is the result of two sign changes, one where a 
— 1 is introduced, and one where a 1 is reintroduced. Descents of the second kind are caused when 
there is no sign change but a descent in tt is present. Descents of the third kind are caused by an 
"incomplete" sign change in e, where a —1 is introduced when reading window notation right to 
left, but not returned to a 1. The equality in (9) follows immediately from these observations. 

To prove (10), observe that the expression 2j contributes to the type- A major index for 

jGDcs(7r) 

all descents that do not involve a sign change. It is then straightforward to see that 

E J = neg(7r, e) + ^ , 

j-aj=l iGDes(7T,e) 

since on the right-hand side, each descent in (vr, e) given by a sign change contributes twice the 
descent position plus the number of consecutive negative signs to the right of the descent position. 
As all negative signs are counted once in this representation, and all descent positions are counted 
once in each of the two summands, the right-hand side of (10) is equal to 2-major^(7r, e)-|-neg(7r, c). 
This completes our proof. □ 

In the original work of Adin, Brenti, and Roichman [1], it was left as an open question to give 
a bijective proof of the equidistribution of the pairs of statistics (ndes, nmaj) and (fdes, fmajor); 
a combinatorial proof leading to an implicit bijection was given by Lai and Petersen in [18]. The 
proofs of Theorems 5.7 and 5.9 together yield this bijective proof. 
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Corollary 5.10 (Adin-Brenti-Roichman) . 

^ ^ ^ndes(7r,e)^nmaj(7r,e) _ ^ ^ ^fdes(7r,e)^fmajor(7r,e) 

(7r,e)GB„ (7r,e)GB„ 

Proof. Our proof arises bijectively from the indexing of integer points in fundamental parallelpipes 
for cone(A^) found in the proofs of Theorems 5.7 and 5.9. To the element (p, e)7r G [JneSn -^2,n'?i" 
we associated the integer point 

n 

p= eo + e^(i) + • • • + e^(j) + ^ 2aj (eo + e^(i) + • • • + e^(j) ) 

jGDes(7r) j=l 

where 

if j^Ncg((p,e)-i), 
1 

2 

Representing p as 



"^■-^ i if iGNeg((p,e)-i^ 



P = H + ^'^(i) + ■ ■ ■ + ^ttO)) + XI ^(^0 + ^'^(i) + ■ ■ ■ + ^ttO)) 

J=l jeDes(7r) 

0j=O 

where Pj € {0, |}, we associated to p the integer point (tt, e) G .Bn where where tt is the same as 
the index on cone(A7r) and e is defined by Qj = 2(3j. This yields an explicit bijection from Bn to 
itself that preserves the pairs of statistics (ndes, nmaj) and (fdes, fmajor). □ 

5.3. g-analogues involving (2A; + 1)". Recall Definition 2.1 which introduced NatDcs(7r,e) and 

natdes(7r, e). Further, for an element (tt, e) G -B„, the naturally ordered major index is 

natmaj(7r, e) := i. 

ieNatDes(7r,e) 

Chow and Gessel [11, Equation (26)] proved the following hyperoctahedral analogue of Theorem 1.3: 
Theorem 5.11 (Chow-Gessel) . 

V c r oneg(e)^natdes(7r,e)„natmaj(7r,e) 

The special case g = 1 is due to Brenti [9, Theorem 3.4]. Chow and Gessel also showed in [11] 
how Theorem 5.11 implies other versions of "g-Eulerian polynomials" of type B involving a flag 
major statistic using the natural order, such as the following. 

Definition 5.12. Use the order —n <■■■<— 1 <!<■■■ <n on [—n,n] \ {0}. We define the 
natural type-A descent set as 

NatDesA(7r,e) := {i e [n - 1] : eiiTi > ej+iTTj+i} 

and the natural type-A major index as 

natmajor^(7r, e) := i. 

iSNatDesA(7r,e) 

The natural flag major index is 

natfmaj(7r, e) := 2 • natmajor^(7r, e) + neg(7r, e) . 
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Theorem 5.13 (Chow-Gessel) . 



EE j.natdes(7r,e) natfmaj(7r,e) 



fc>0 



While Theorems 5.11 and 5.13 are equivalent via a change of variables, the geometric perspective 
illustrates how these two theorems arise as specializations of two distinct multivariate generating 
function identities. 

5.4. Multivariate identities. For the type-i? extension of Theorem 4.1 generalizing Theorem 5.11, 
we introduce the variables z±j to keep track of the positive/negative jth component of a lattice 
point, respectively, and the variables w±j to indicate whether the _7th coordinate of our point is 
nonnegative or negative. 

Theorem 5.14. 

n 

Yi n + + ^-j ^-ji^izj) ^0 = 

fe>0 j=l ^ 



ET-r T-r _i jrGNatDes(7r,e) 
[[Wj ii ^-jW-J^ -• 

j=0 

Note that Theorem 4.1 is the special case w-i = ■ ■ ■ = w-n = 0, lui = • • • = = 1 of 
Theorem 5.14. 

Proof. Triangulate the n-cube [—1, 1]" by first disjointly decomposing [—1, 1]" into the cubes 

y y) < EjXj <1 if Ej = —1 

where there is one cube for every e G {±1}". Each of these cubes then gets triangulated into the 
disjoint simplices 

^ •= /a; G R" • ^ ~ ^"^^(n) < en-ia:7r(n-i) < • • • < eia;7r(i) - ^ 1 c □ 
'''''' 1 ' ej+ia;^(j+i) < ejx^(j) if j G NatDes(7r,e) J - ^ 

(recall that we set eo = 0). A lattice point m G cone(ne) gets encoded by the monomial 

zr n --.-5 n-.-r - 
ei=-i ej=i 

Note that for each j, exactly one of Zj and z-j appears in this formula. Since cone(A7r,e) C cone(ne) 
(and so for every j, the sign of the jth coordinate of any point in cone(A7r,e) is the same), we can 
adapt Lemma 3.1 to see that 



){ZQ,Z±i,...,Z±n,W±i,...,W±n) = 

''ei7r(l)'^e27r(2) '^'^j'^ij) 



n n n - 



jeNatDes(7r,e) 



j=0 
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Again note that for each j, either zj or z-j appears in this formula, but not both. On the other 
hand, 

'7conc([-l,l]")(^0, Z±i, . . . , Z±n,W±i, . . . ,W±n) = 
n 

n {^-j ^-j + ^-f~^^ + • • • + zZ] + Wj + Wj Zj +Wjz] + --- + Wj z^ , 

k>Oj=l 

and the disjoint triangulation gives 

O"cone([-l,l]")(^0, Z±l,..., Z±n, W±l, W±n) = ^ O-cone(A^) (^0, Z±l, . . . , Z±n, W±l, W±n) ■ 

(7r,e)GBn 

□ 

Proof of Theorem 5.11. Setting t := Zq, q := Z\ = ■ ■ ■ = z^ = zZ\ = ■ ■ ■ = zZ^, s := w^i = ■ ■ ■ = 
w-n, and lui = • • • = ■u;^ = 1 in Theorem 5.14 gives 

E n^^ n V 

E„neg(e)j.natdcs(7r,e) natmaj(7r,e)+neg(7r,e) 
^ (7r,e)eB„ ^ ^ y 

n -=0 (1 - tq^) 

Replacing s by | yields Theorem 5.11. □ 

Our multivariate generalization of Theorem 5.13 is the following. Recall from Definition 5.8 the 
notation ch(e) for the number of sign changes in e and the notation aj to keep track of where sign 
changes occur. 

Theorem 5.15. 



En [2^+1]^ 

A;>0 j=l 

^rch(e)/2l TT^,,^,,...^,, TT ^„^2 ^2 ... ^2 

j-a^j=^ jeDes(7r) 
V- •'j^ 

/ y n 

Proof. We begin by triangulating [0, 2]" into the disjoint simplices 

1^ 2;^(j+i) < if J e Des(7r) 

(one for each vr G S'n). As before, the strict inequalities determined by the descent set of tt ensures 
that this triangulation is disjoint. 

Unlike the cones produced by coning over the simplices in our triangulation of [0, 1]", the cones 
arising from this triangulation of [0, 2]" are not unimodular. By Lemma 3.2, the integer-point 
transform of cone(A7r) can be expressed as a rational function with denominator 



n(^~^o4(i)4^''^---^' 



4(2) • • • 4(j) ; > 
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i.e., where the displayed exponent vectors are the ray generators for this cone. As the determinant 
of the matrix formed by the ray generators of cone(A^) is 2", there are 2^* integer points in 
the fundamental parallelepiped of cone(A7r). It is a straightforward observation that there are 
2" such integer points formed by taking linear combinations of the ray generators for the cone 
with coefficients and ^. We will use the following notation to denote such an integer point; for 
a, e {0,i}, 

n 

p = aoeo + aj (eo + 2e^(i) + 2e^(2) + • • • + 2e^(j)) . 
i=i 

To determine the numerator of o'cone(A^)(-20) • • • j Zn), we must shift some of these interior points off 
of the boundary of cone(A7r). It is a simple exercise to show that the point p must be shifted by 
the vector 

^ eo + 2e7r(i) + 2e^(^2) H ^ 2e^(j) . 

j:aj=0 

Thus, 

O'cone(A„)(^0, Zi, . . . , Zn) = 

nil 2 
^0^7r(l)^7r(2) ■■■^7r(j) 

j€Des(7r) 

Qj=0 

n 

n (i-^o 4(1)4(2) 

i=o 

Note that G [0, n] : = ^} is always an integer, as this is the coefficient of eo in the resulting 

integer point. 

We now associate to the element (tt, e) G -B„ the integer point 

" a"- 

J=l jeDes(7r) 

n 

in cone(A^), where ao = or ^ as needed to make ao + X] integer. Through this association, 

the set of sign vectors {±1}" parametrizes the integer points in the fTindamental parallelepiped for 
cone(A^). This parametrization is bijective, and the coefficient of eo in the first two terms of the 
sum above is equal to both [ch(e)/2] and € [0, n] : = ^}. Thus, our proof is complete 

following the observation that 

f^cone([0,2]")(^0, • • • , ^^n) = X '^cone(A„) (^0, • • • , ^n) • □ 

■KeSn 

Proof of Theorem 5.13. The pairs of statistics (natdes, natfmaj) and (des, fmajor) are equidistributed 
in Bn, this is a consequence of the bijection mapping every permutation (tt, e) G B„ to the permu- 
tation where the 7r(j) for ej = —1 have been reversed in order. Hence, we may prove our theorem 
using the second pair of statistics. 

Setting t = zq and g = zi = • • • = z„ in Theorem 5.15 yields our desired form on the left- 
hand side, while the denominator of the right-hand side uniformly becomes 11^=0 ~ ^(f'')- Each 



i#{ie[0,n]:a, = i} 



n 



Z'jv{l)ZTr[2) ■ ■ ■ 



3:aj-- 
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element (tt, e) G Bn contributes a summand of 

^[ch(e)/2l -Q Yl tq^^ 

i-«5=l ieDes(7r) 

to the numerator of the right-hand side. Hence, we need to prove 

(11) des(7r, e) = [ch(e)/2l + #{j € Des(7r) : a'j = 0} , 
and 

(12) fmajor(7r,e)= ^ j+ 2j . 

j:a^j=l jeDes(7r) 
a^.=0 

As explained in the proof of Theorem 5.5, a descent in position j of (vr, e) can arise for one of 
three reasons: either ej = 1 while e^+i = —1, or ej = tj+i and j € Dcs(7r), or j = 1 and ei = — 1. 
When viewing (tt, e) in window notation and reading right to left, every descent of the first kind is 
the result of two sign changes, one where a —1 is introduced, and one where a 1 is reintroduced. 
Descents of the second kind are caused when there is no sign change but a descent in vr is present. 
Descents of the third kind are caused by an "incomplete" sign change in e, where a —1 is introduced 
when reading window notation right to left, but not returned to a 1. The equality in (11) follows 
immediately from these observations. 

To prove (12), observe that ^ 2j contributes to the type-^ major index for all descents 

jeDes(7r) 

that do not involve a sign change. It is then straightforward to see that 

^ i = neg(7r,e)+ ^ 2 j , 

j:a^=l jeDcs(7r,e) 

since on the right-hand side, each descent in (tt, e) given by a sign change contributes twice the 
descent position plus the number of consecutive negative signs to the right of the descent position. 

As all negative signs are counted once in this representation, and all descent positions arc counted 
once in each of the two summands, the right-hand side of (12) is equal to 2 • maj^(7r, e) + neg(7r, e). 
This completes our proof. □ 



6. Type D 

6.1. (/-analogues involving {k + 1)". Moving on to the groups Dn, we recall the following defi- 
nition of type-I? negative statistics [6] . 

Definition 6.1. Using the order —1 < • • • < —n < 1 < • • • < n on [—n,n] \ {0}, for an element 
(tt, e) G Dn, we define Neg(7r, e), neg(7r, e), Des^(7r, e), and des^(7r, e) as for the group Bn- We 
define the type-D negative descent multiset as 

DNDes(7r, e) := Bcsa{tt, e) [j{n{i) - 1 : q = -1} \ {0} 

= BesAin, e) \J{j - 1 : j e Neg((7r, e)-') \ {1}} . 

The type-D negative descent statistic is 

dndes(7r, e) := #DNDes(7r, e) . 
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The type-D negative major index is 

dnmaj(7r, e) := i. 

jeDNDcs(7r,e) 

Again in this case, Biagioli [6] defined dnmaj and dndes using the natural order, but the following 
theorem holds for either definition. 



Theorem 6.2 (Biagioli). 



j.dndes(7r,e)„dnmaj(7r,e) 

y[k+ 1]- 1'^ = 



6.2. Multivariate identities. Our multivariate generalization of Theorem 6.2 is as follows. Let 
C l2^n denote the elements (p, e) G l2,n satisfying €±€2 • • • e„ = 1. It is straightforward to show 

that 

Dn= [j /2>- 

Theorem 6.3. 

n 

fc>oi=i 

^ jGDes(7r) jeNeg((p,c)-i)\{l} 



E 



^^s^ M^iln (1 _ _ . . . ^^^^^) JJ (1 _ ^2^2^^^ . . . 4^ 

Proof. We begin with the triangulation of cone([0, 1]") into the set of cones {cone(A7r) : tt G 5^} 
found in the proof of Theorem 4.1. For cone(A7r) we use the non-unimodular ray generators 

eo, 2(eo + e^(i)), . . . , 2(eo + e^(i) H h e7r(„_i)), cq + e^(i) H h e^(„) . 

There are 2"'~-'^ integer points in the fundamental parallelepiped for cone(A^) using these ray 
generators. Each such point can be expressed as a linear combination of the middle n — 1 generators 
with coefficients aj G {0, |}, plus a sum of shifting vectors for those integer points that need to 
be shifted away from the boundary of the cone. As in our proof of Theorem 5.7, we will shift all 

the points in the fundamental parallelepiped by the minimal ray generators eo + ^^{i) H + ^wij) 

when j G Des(7r). This can again be justified by noting that the fundamental parallelepiped for 
cone(A7r) with our non-unimodular ray generators is lattice-isomorphic to a half-open cube of the 
form [0,1) X [0,2)"-i x [0,1). 

Thus, every point p in the (shifted) fundamental parallelepiped for cone(A7r) can be uniquely 
expressed as 

n-l 

p= E eo + e^(i) H h e^Q) + ^ 2aj (eo + e^(i) H h e^(j)) 

jeDes(7r) j=l 

with aj G {0, ^}. Associate to the point p the element (p, e)7r G D„, where aj = 5 if and only 
if J -|- 1 G Neg((/9, e)~"^) = Neg([(/9, e)7r]~^). This correspondence creates a bijcction between the 
elements of Dn and the (appropriately shifted) integer points in the fundamental parallelepipeds for 
the cones over the A^^; in this case the bijection is less obvious than in the case for Our choice of 
(p, e)7r associated to p is unique because the condition aj = 5 if and only if j -|- 1 G Neg((p, e)""^) = 
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Neg([(p, e)7r]~^) determines the signs placed on the letters {2, 3, ... , n} when (p, e)7r is written in 
window notation. Hence, ei is determined from these n— 1 signs and the fact that ei • • • e„ = 1. Note 
as before that this bijection encodes as the integer points in the fundamental parallelepiped 
for cone(Ai(j). Thus 

IJ ^0^7r(l)2^7r(2) • • • Zj^ij) Yl ^QZ'k{1)Z'k{2) ' ' ' ^ttO-I) 

jeDes(7r) jeNeg((|0,e)-i) 

<''cone(A^)l^O, • • • ,2;nj = 2^ 



Meiin (1 _ - zoz,(i) • • • n (i - ^0^(1) • • • ^. 



n-1 

2 



This completes our proof, since from our triangulation it follows that 

''^cone([0,l]") ^cone(A7r ){Z0, ■ ■ ■ , Zn) ■ □ 

Proof of Theorem 6.2. Setting t = zq and q = zi = ■ ■ ■ = Zn m Theorem 6.3 yields our desired 
form on the left-hand side of our identity, while the denominator of the right-hand side uniformly 
becomes 

n-1 

(l-t)(l-tg")n(l-*V^)- 

Each element (p, e)n G /2,n7i" contributes to the numerator on the right-hand side of our identity 
a summand of 

n n 

jeDes(7r) ieNeg([(p,e)7r]-i) 

Because Des(7r) = DesA((/9, e)7r), it follows that 

n n tgJ-l = idndes(p(7r,e))^dnmaj(p(7r,e))^ 

jeDes(7r) jeNeg([(p,e)7r]-i) 

hence our proof is complete. □ 

6.3. g-analogues involving (2A;+ 1)". We derive in this subsection a type-D analogue of Theorem 
5.11, in the sense that we'd like to have the same left-hand side and the same denominator on the 
right-hand side. A natural starting place in this endeavor is to seek a (/-analogue of Brenti's identity 
(7). We should mention that in [9], Brenti stipulated that he proved an "s-analogue" of (7) (where 
s plays the same role as in Theorem 5.11), while Biagioli [6, Theorem 3.4] and Biagioli-Zeng [7, 
Theorem 6.3] computed different type-D analogues. 

In type D, the (possible) descent positions and 1 are special; thus we define 

N2{e) :=#{e, : j > 2, = -1} , 
dnatdes2(7r,e) := # (DNatDes(7r, e) n [2,n]) , 

maj2(7r,e):= ^ j. 

jgDNatDcs(7r,e)n[2,n] 

Recalling that Bn{k) denotes the v}^ Bernoulli polynomial and referencing (7), since 

k 

Bn{k + l)-Bn{Q) = nY,f-\ 

3=0 
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we can simplify 

Y: 2^-^ {Bnik + 1) - B^m = ^ ^0- + 1)"- = ll'il'-i — • 

Via the natural (n 2"~^)-to-one map D„ Sn-i given by (7r,e) i-> (e27r(2), e37r(3), . . . , eH7r(n)), 
where we interpret (e27r(2), e37r(3), . . . , e„7r(n)) as an element of Sn-i, we obtain 



dnatdes2(7r,e) 

7re5n-i (7r,e)eD„ 



This gives the following reformulation of (7), which — as we will see — arises naturally from our 
approach: 



j.dnatdi 

(13) ^(2A: + lf 

For a proposition P, we use the notation 



^dnatdes(7r,e) _j_ ^l+dnatdes2(7r,e) 

n+1 



[P\ : - 



1 if P is true, 
if P is false. 



Our g-version of (13) is the following. 
Theorem 6.4. 

Y{[k+iu+s[k],rt' 

k>0 

Z](7re)eD S-'^2(e)^dnatdes2(7r,e)^maj2(7r,e) ^^^^yei=-l] ^^^yieDNatDes(7r,e)] _|_ ^^^^[ei=l] ^^^^[06DNatDes(7r,e)]^ 

^ n-=o(i-i?^) 

y^ ^ ^ ^ gAr2(e)^dnatdes2(7r,e)gmaj2(7r,e) ^^^^yOor leDNatDes(7r,e)] _j_ g^^^g^[Oand leDNatDes(7r,e)]^ 

n-=o(i-*9^) ■ 

Brenti's identity (13) follows upon setting s = q = 1 and noticing that 

^dnatdes2(7r,e) ^^[Oor leDNatDes(7r,e)] _j_ ^ _ ^[OandleDNatDes(7r,e)]^ ^dnatdes(7r,e) _j_ ^l+dnatdes2(7r,e) 

As far as we can tell, no such simplification happens when we use a different specialization for s 
and q. 

Proof of Theorem 6.4- For a given (tt, e) € Z)„, let e := (— ei, €2, . . . , e„). Then 

(7r,e)gD„ 

is the hyperoctahedral group, and so by Theorem 5.11, 

('„iV(e)/natdes(7r,e) „maj(7r,e) i „7V(e)j.natdes(7r,€) „maj(7r,e)'\ 

E {[k + 1], + . [k],r = ^^-'-^^""^ — ' — ^ . 
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In terms of these descents, (tt, e) and (tt, e) differ only (possibly) in positions and 1, and so a 

summand in the numerator looks like 

(14) 

giV(£)^natdes(7r,e)gmaj(7r,e) _|_ gAr(e)^natdes(7r,e)^maj(7r,e) 

_ ^iV2(e)^dnatdes2(7r,e)^maj2(7r,e) f |-^^-j[0eDNatDes(7r,€)] ^^^-j [leDNatDes(7r,e)] _|_ ^^^y0eNatDes(7r,e)] |-^^yiGNatDes(7r 



The role of s in the large parenthesis comes from the fact that ei = — 1 if and only if € NatDes(7r, e). 
Furthermore, since 

1 G NatDes(7r, e) <S=^ eiTTi > e27r2 1 G DNatDes(7r, e) 

and 

1 G NatDes(7r, e) <;=^ -citti > €21^2 G DNatDes(7r, e) , 

we can rewrite (14) as 

^Ar(e)^natdes(7r,e)^maj(7r,e) _|_ gAr(e)^natdes(7r,e)^maj(7r,e) 

^ gAf2(e)^dnatdes2(7r,e)^maj2(7r,e) ^(^g^^[ei=-l] |-^^yiGDNatDes(7r,e)] _|_ ^^^-)[ei=l] (^^^yOGDNatDes(7r,e)]^ _ g 

7. Wreath products 

7.1. g-analogues involving (A; + l)". We now consider statistics for the groups I Sn, these 
extend and subsume the statistics introduced for the typc-i? case. We begin by defining typc-yl 
statistics for elements of I Sn, [3, 4]. Throughout, we use the total order from Definition 2.3 on 
the elements of {uj^~^ , u!^~'^ , . . . ,uj^} x [n], i.e., j'^^ < k'^'^ if Cj > cj- or if both cj = and j < k 
hold. 

Definition 7.1. For an element (tt, e) El>rlSn, we define the negative set of (tt, e) to be 

Neg(7r, e) := {i G [n] : 7^ a;° = 1} , 
and we define neg(7r, e) := #Neg(7r, e). Writing Sj = oj'^^ , we define the color sum statistic to be 

col(7r,e) := ^ 

The type-A descent set is defined to be 



Ci 
ie\n\ 



DesA(7r, e) := {i G [n - 1] : Trf > 
and the type-A descent statistic is 

desA(7r, e) := #DesA(7r, e) . 
We next define negative statistics for wreath products, following [3]. 
Definition 7.2. For an element (7r,e) in I Sn, we define the negative inverse multiset as 

NNeg(7r,e) := {i"' : i G [n]}. 
We define the negative descent multiset as 

NDes(7r, e) := DesA(7r, e) (J NNeg((7r, e)"^) . 
The negative descent statistic is 

ndes(7r, e) := #NDes(7r, e) . 
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The negative major index is 

nmaj(7r,e) := ^ i. 

iGNDes(7r,e) 

Observe that NNeg((7r, e)^^) contains exactly (r — c^-i(j)) mod r copies of each i £ [n]. The 
negative statistics for wreath products lead to the following theorem due to Bagno [3]. 



Theorem 7.3 (Bagno). 



EE j.ndes(7r,e) „nmaj(7r,e) 
,J (l-i)n"=i(l-i^?^^) 



fe>0 

There are also flag statistics for wreath products, due to Bagno and Biagioli [4] . 

Definition 7.4 (Bagno-Biagioli) . For an element (tt, e) in Zj-lSn, we define the flag descent statistic 
as 

fdes(7r, e) := r • des^(7r, e) + ci , 
where as usual ei = u)^'^ . The flag major index is 

fmajor(7r, e) := r ■ major^(7r, e) + col(7r, e) . 
The flag statistics for wreath products lead to the following theorem. 
Theorem 7.5 (Bagno-Biagioli). 



EE j.fdes(7r,e) fmajor(7r,e) 



Bagno and Biagioli prove in [4] a multivariate theorem of this type for a family of normal 
subgroups of Zr I Sn- Their techniques involve studying colored-descent representations of these 
subgroups, which are representations of the groups on the associated coinvariant algebra. 



7.2. Multivariate identities. We state in this subsection new multivariate generalizations of 
Theorems 7.3 and 7.5. When our theorems arise from straightforward modifications of our proofs 
in type B, we will sketch or omit their proofs. 

Definition 7.6. Define the subset of increasing elements of I Sn, denoted Ir,n, to be those 
elements satisfying des^(p, e) = 0, i.e., Ir,n contains all permutations (p, e) such that p{jY^ < 
p{j + lp+K 

It is straightforward that every element of "LylSn can be represented uniquely as (p, e)7r = (p7r, e-,^) 
for some n E Sn and (p, e) G Ir,n) where pir is the standard composition of permutations and 

^TT = (e7r(l)) e^(2)) ■ ■ ■ ^77{n))- Thus, 

7re5„ 

Observe that Neg([(p, e)7r]~^) = Neg((p, e)~^). Our multivariate version of Theorem 7.3 is the 
following. 
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Theorem 7.7. 

n 

fe>0 j=l 

jGDes(7r) iGNcg({p,e)-l) 



E E 



(^'^)^^- (1 - n (i - ^5<(i) • • • 

i=i 

The following definition is needed to state our multivariate version of Theorem 7.5. 

Definition 7.8. For an element e = (w'^j)^L;^ G {1, w-^, w^, . . . , w*'^-*^}" with e„+i := 1 = a;°, for 
j G [n] define 

a^- := {cj — Cj+i) mod r , 
which we call the jth sign change for e. Define ch(s) := a^- to be the toto/ sign change in e. 

Theorem 7.9. 



n 

En(i+^^+4+---+^')^o 



A;>0 j=l 

n 



n ^o4(i) 4(2) • • • n V(i) V(2) ■ ■ ■ 

jeDes(7r) j=l 

E 



aj=0 



(.,e)e^.i5„ (1 _ n (l - <(1)<(2) • • • 

The proofs of Theorems 7.7 and 7.9 together yield a bijective proof of the equidistribution of 
the pairs of statistics (ndes, nmaj) and (fdes, fmajor) for I Sn- As far as we know, this bijection 
is new. The proof of this bijection follows immediately from the proofs of Theorems 7.7 and 7.9, 
whose proofs are easily obtained by modifying the proofs of Theorems 5.7 and 5.9. 

Corollary 7.10. 

^ ^ ^ndes(7r,€)^nmaj(7r,e) ^ ^ ^fdes(7r,e)^fmajor(7r,e) 

7.3. g-analogues involving arbitrary (r/c + 1)". In [12, Theorem 9], Chow and Mansour provide 
an Euler-Mahonian distribution for wreath products using Steingn'msson's descent statistics and 
a new flag major index. We state in this section a similar Euler-Mahonian distribution for the 
descent statistic and flag major index given in Definitions 2.3 and 7.4. By combining Theorem 7.11 
below and [12, Theorem 9], we see that the pairs 

(Steingrfmsson's descent statistic, Chow-Mansour's flag major index) 

and 

(des, fmajor) 

are equidistributed over "LylSn- 
Theorem 7.11. 

j.des(7r,e) fmajor(7r,e) 

Ylrk + 1]- = ^^"-^^"fcf ' . 
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We obtain in Theorem 7.13 below a multivariate generalization of our bivariate identity. At the 
end of [12], Chow and Mansour indicate that having a Hilbert-series interpretation of [12, Theorem 
9] is an open problem. The proof of Theorem 7.13 provides such an interpretation, after taking 
into account Remark 4.3. 

Remark 7.12. In [8, Equation (8.1)], Biagioli and Zeng obtain a gencrahzation of the Chow-Gessel 
Theorem 5.11 for wreath products. The authors do not at this time know of a way to obtain this 
identity using polyhedral geometry. 

7.4. Multivariate identities. 

Theorem 7.13. 



fc>oi=i 

n 

rch(6)/rl-|-r a' yr 

j=l jeDes(7r) 
'^J^ 

Proof. As in our previous proofs, this proof proceeds in two stages. We first triangulate the cube 
[0, r]"^ into a disjoint union of simplices, then set up an indexing system for the integer points in the 
fundamental parallelepipeds for the cones over these simplices. Second, we bijectively associate the 
elements of I Sn with these integer points in a way that allows us to recover, in our subsequent 
proof of Theorem 7.11, the descent and flag major statistics from these integer points. 
We begin by triangulating [0,r]" into the disjoint simplices 

A ■=lxeW- ° - '"""("^ - - ■ ■ ■ - - 

^ ■ \ ■ a;^(j+i) < x^f^j) if j € Des(7r) 

(one for each tt G Sn)- As before, the strict inequalities determined by the descent set of tt ensures 
that this triangulation is disjoint. 

Unlike the cones produced by coning over the simplices in our triangulation of [0, 1]", the cones 
arising from this triangulation of [0, rj" are not unimodular. By Lemma 3.2, the integer-point 
transform of cone(A^) can be expressed as a rational function where the denominator has the form 

n 

n (i-^o4(i)4(2)---4o)) ' 

i=o 

i.e., where the displayed exponent vectors are the ray generators for this cone. As the determinant 
of the matrix formed by the ray generators of cone(A^) is r", there are r" integer points in the 
fundamental parallelepiped of conc(A7r). It is a straightforward observation that there are r"^ 
such integer points formed by taking linear combinations of the ray generators for the cone with 
coefficients from the set {O, ^, . . . , ^-7^}- We will use the following notation to denote such an 
integer point; for G {O, ^, . . . , } where j = 0, . . . ,n, 

n 

p = aoeo + ^ aj (eo + 7-67^(1) + re^(2) H h r-e7r(j)) . 

i=i 

To determine the numerator of crcone(A,r)(^0) • • • , ^n), we must shift some of these interior points off 
of the boundary of cone(A7r). It is a simple exercise to show that the point p must be shifted by 
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the vector 

eo + re^(i) + re^2) + ■■■ + re^{j) ■ 

jGDes(7r) 
aj=0 

Thus, 

'^cone(A,r)(^0) Zi, . . . , Zn) = 

2.^=0 TT raj raj raj yr r r . . r 



'7r(j) 



j=l jeDes(7r) 
ai=0 



n (i - 20 <(i)<(2) •••<(,)) 

Note that Yl^=o is always an integer, as this is the coefficient of Cq in the resulting integer point. 
We now associate to the element (tt, e) G ? Sn the integer point 

" a' 

aoeo + Y^ieo + re^{i) + re^{2) + ■■■ + re^{j)) + Yl + re^{i) + re^{2) + ■■■ + re^{j)) 

J=l ieDes(7r) 

a'j=0 

n e 

in cone(A7r), where ctg takes the value needed to make uq + an integer. Through this 

association, the set of sign vectors {1,^,^"^ , . . . ,uj^~^}"' parametrizes the integer points in the 
fundamental parallelepiped for cone(A^). This parametrization is bijective, and the coefficient of 
eo in the first two terms of the sum above is equal to both [ch(€)/r] and J2]=o Thus, our proof 
is complete following the observation that 

^cone([0,r]")(20, • • • , ^n) = ^ O-cone(A^) (^0, • • • , ^n) • D 

7re5„ 

Rather than include a detailed proof of how Theorem 7.11 follows from Theorem 7.13, we will 
instead sketch the proof and then include an illustrative example. 

Proof of Theorem 7.11. Setting t = zq and g = zi = • • • = z„ in Theorem 7.13 yields our de- 
sired form on the left-hand side, while the denominator of the right-hand side uniformly becomes 
nj=o ~ t'f'')- Each element (vr, e) ^2,rl Sn contributes a summand of 

j-.a^-^O jGDes(7r) 

to the numerator of the right-hand side. Hence, we need to prove 

(15) des(7r, e) = rch(e)/r] + #{i G Des(7r) : a) = 0} , 
and 

(16) fmajor(7r, e) = XI ^^3^ + XI ' 

3-a-i¥'0 j6Des(7r) 

As explained in the proof of Theorem 5.5 in the case of type B, a descent in position j of (tt, e) can 
arise for one of three reasons: either Cj > Cj+i, or Cj = Cj+i and j G Des(7r), or _7 = 1 and ci 7^ 0. 
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When viewing (tt, e) in window notation and reading right to left, every descent of the first kind 
is the result of successive sign changes that sum to a value greater than or equal to r. Descents of 
the second kind are caused when there is no sign change but a descent in vr is present. Descents 
of the third kind are caused by an "incomplete" sign change in e, i.e., where non-zero colors are 
introduced when reading window notation right to left, but these colors are not returned to a 
before reaching the first position on the left. The equality in (15) follows immediately from these 
observations. 

To prove (16), observe that rj accounts in the type- A major index for the contribution 

jgDcs(7r) 
a5=0 

given by those descents that do not involve a sign change; it then only remains to consider the 
contribution of those descents arising from a sign change. Therefore, what remains to be shown is 
that 

(17) a'jj = col{Tr,e)+ ^ rj. 

3-a^j¥=0 jeDcs(7r,e) 

We will sketch a proof of this and include an example illustrating it. 

In order for a descent to arise as a sign change, it must have been that reading right to left in 
window notation for {TT,e), the sum of the sign changes to that point must have exceeded r — 1, 
hence returning the color for the descent position to a value less than the one preceding it on the 
right. Therefore, if the descent occurred in position d, the term 

Y. "V 

will have contributed rd for that descent plus some additional amount. This additional amount 
then becomes distributed between the term col(7r, e) for the colors q where i > d, and the sum rd' 
for the next descent position d' to the left of d. Hence, (17) holds, and our proof is complete. □ 

Example 7.14. To illustrate the final step in the proof above, we will show that for the color vector 
c = 412301131 arising in a possible element of Z5 I Sg, (17) holds. Note that there are four type-^ 
descents in c caused by sign changes, occuring in positions 7, 5, 3, and 2. The sign change vector 
for c is a = 344340321, where the right-hand 1 accounts for cg — cio where cio is by definition 0. 
We then have that 

9ag + Sag + 7ar = 7- 5-|-C9-|-C8 + 7. 

Hence, the descent in position 7 contributes 7 • 5 -|- cg -|- cg to the left-hand side of (17), plus 7 more 
which is counted in 

7 + 5a5 = 5 • 5 + C7 + C6 . 
This yields the contribution of the 5th descent position plus colors cq and cy to the left-hand side 
of (17). Next we obtain that 

4a4 -|- 3a3 = 3- 5-|-C5-|-C4-|-6, 
contributing the 3rd descent position plus colors C5 and C4 to (17) with 6 remaining. This leftover 
6 then contributes to 

6 + 2a2 = 2 • 5 + C3 + 2 , 
where the final 2 in the expression above contributes to 

2 -I- lOl = C2 + Cl . 

Thus, the sum of juj over j is equal to the statistic col(7r, e) plus the terms rj in the type-A major 
index statistic arising from sign changes rather than descents in tt. 
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